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Abstract 

The  primary  source  of  aircraft  noise  is  the  fan  noise  from  the  engines;  natural  approaches 
to  reducing  this  noise  involve  acoustic  shape  optimization  of  the  inlet  and  impedance  opti¬ 
mization  of  the  liner.  This  project  will  use  optimal  control  to  systematically  determine  the 
inlet  shape  and  the  linear  material  impedance  factor  that  minimize  the  fan  noise. 

A  novel  feature  of  this  approach  is  that  we  automatically  incorporate  uncertainty  and  data 
measurement  errors.  Specifically  we  assume  that  the  acoustic  wave  number  is  a  random 
variable/field  instead  of  a  constant.  This  means  that  the  computed  answers  are  valid,  not 
merely  for  a  single  configuration,  but  for  a  wide  range.  Our  numerical  results  show  signifi¬ 
cant  noise  reduction  with  the  optimal  impedance  factor. 

Since  the  wave  number  is  random,  the  underlying  partial  differential  equation-Helmholtz 
equation  in  our  case,  is  a  stochastic  partial  differential  equation.  In  this  project,  we  have 
constructed  efficient  Monte  Carlo  methods  as  well  as  stochastic  finite  element  methods  to 
solve  stochastic  partial  differential  equations.  Rigorous  error  estimates  are  obtained  and 
numerical  simulations  are  conducted  to  support  the  error  analysis. 

Formulation  of  the  optimal  control  problem 

In  this  research,  we  treat  liner  impedance  optimization  as  an  optimal  control  and  parameter 
estimation  problem.  The  parameter  is  the  acoustic  impedance  factor  of  the  acoustic  liner. 
We  define  a  cost  function  that  reflects  the  amount  of  noise  radiated  from  the  engine  inlet. 
The  parameter  estimation  problem  then  is  to  seek  the  parameter  that  minimizes  the  cost 
function.  The  geometry  of  the  domain  in  which  the  control  problem  is  posed  has  the  generic 
shape  represented  in  Figure  1.  The  modal  composition  of  the  noise  source  is  supposed  to 
be  known  on  the  source  plane  Ti.  The  nacelle  boundary  is  made  up  of  two  parts,  the  first 
part  being  the  interior  boundary  T2  to  which  some  acoustic  liner  material  is  attached,  and 
the  second  part  being  T3  that  constitutes  the  rest  of  boundary  of  the  nacelle  geometry.  The 
boundary  Tzt  is  assumed  to  be  sufficiently  far  from  the  noise  source  so  that  the  Sommerfeld 
radiation  boundary  condition  holds.  The  nacelle  symmetry  axis  is  denoted  by  r5. 

We  assume  that  the  mean  flow  is  zero.  Then  the  acoustic  pressure  u  satisfies  the  Helmholtz 
equation 


A u  +  k2u  =  0  on  9 
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Figure  1.  The  computational  domain 


subject  to  the  following  boundary  conditions  on  the  boundary  dQ  of  Q: 


(2) 


u  |n  =  g, 

+  =  °> 

S  lr3  =  0, 

ffi  +  iku )  lr4  =  0, 

L  _  n 


Here  k  is  the  wavenumber  and  £  the  impedance  factor,  which  is  complex,  and  whose  real  part 
is  resistance  and  the  imaginary  part  is  reactance.  Both  the  dependent  and  the  independent 
variables  in  the  above  equations  are  supposed  to  be  properly  non-dimensionalized. 


We  first  assume  that  k  is  a  constant.  Then  optimization  problem  consists  in  finding  the 
parameter  £  so  that  the  least  amount  of  noise  propagates  to  the  far  held,  while  respecting 
some  constraints  on  the  boundary  shape.  More  specifically,  we  want  to  find  £  so  that  the 
cost  functional 

(3)  J(£,  u)  —  a  [  u2dQ  +  P  [  \Vu\2dfl  +  A|£  -  £0|2 

Jn  Jn 

is  minimized,  where  a,  (5  and  A  are  penalty  parameters,  a  >  0,  f3  >  0,  A  >  0,  and  £o  is  a 
given  complex  number. 

Because  of  the  varying  operating  conditions  and  errors  in  measurement,  k  is  not  likely  a 
constant,  instead  it  may  be  a  random  held.  For  simplicity  we  assume  that  it  is  a  random 
variable  with  small  variance: 

k  =  k0  T  6T] 

where  ko  is  a  constant  and  //  is  a  random  variable  with  zero  expectation,  i.e.,  Etj  =  0.  In 
this  case  we  construct  the  new  cost  function  as 

The  optimal  control  problem  becomes  hnding  %*  such  that 

J{  X*)  =  min  J(x). 

Clearly  such  an  approach  makes  the  control  models  more  complex,  but  much  more  hexible, 
realistic  and  practical.  Because  the  effects  of  data  uncertainty  are  built  into  the  model,  we 
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expect  to  see  that  the  optimal  controls  will  be  much  less  sensitive  to  changes  in  the  model 
parameters.  The  result  is  much  greater  robustness  for  the  optimal  controls. 

Fast  Monte  Carlo  simulation  for  the  optimal  control  problem 

The  main  challenge  of  solving  the  optimal  control  problem  under  uncertainty  is  the  need 
of  solving  large  number  the  state  equations  when  using  Monte  Carlo  method  to  evaluate 
the  cost  function.  We  have  developed  a  fast  numerical  algorithm  to  obtain  the  numerical 
solution  for  the  optimal  control  problem  under  uncertainty  ([1,  2])  .  The  main  idea  is  to 
reduce  the  variance  by  sampling  the  residual  of  the  Taylor  expansion  for  cost  function  in 
terms  of  the  random  wave  number. 

We  have  performed  a  variety  of  numerical  experiments  to  show  that  the  choices  of  optimal 
parameters  result  in  reducing  the  noise  level  significantly.  Here  we  only  present  one  example 
to  illustrate  the  numerical  procedure.  We  choose  the  wave  number  Ay,  =  2tt,  e  =  0.2  and 
r]  a  random  variable  with  uniform  distribution  U(— 1,1).  The  source  function  is  defined 
by  g(y )  =  20  +  exp(3y)sin(10Try).  To  evaluate  the  cost  function,  we  combine  the  sensitive 
derivative  Monte  Carlo  (SDMC)  developed  in  [1]  with  the  stratified  Monte  Carlo  method  to 
accelerate  the  evaluation  of  the  cost  function  (see  [2]  for  detailed  implementation).  Figure 
1  is  the  contour  map  of  the  attitude  of  the  acoustic  pressure  before  and  after  control.  It 
indicates  significant  noise  reduction  with  optimal  impedance  factor. 

The  contour  maps  of  the  amplitude  of  the  acoustic  pressure  u  are  displayed  in  Figure  3  for 

£  =  oo  and  optimal  impedance  factor  £  =  £*. 


Figure  2.  The  contour  maps  of  |£’(u(oo))|  and  |i£(w(£*))|:  a  =  1,  (3  =  1 


Numerical  solutions  for  the  underlying  stochastic  partial  differential  equations 
(SPDES) 

Because  of  the  randomness  of  k,  the  Helmholtz  equation  is  a  stochastic  partial  differential 
equation.  The  Monte  Carlo  method  is  a  generally  applicable  solution  method  for  stochastic 
problems.  However,  it  is  a  well-known  fact  that  as  the  accuracy  requirement  is  increased, 
the  number  of  realizations  to  be  generated  and  deterministic  problems  to  be  solved  grows 
far  too  rapidly.  In  [3],  We  constructed  a  finite  element  method  based  on  polynomial  chaos 
expansion  (PCE)  and  obtained  the  following  rigorous  error  estimate. 
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Theorem  1.  Assume  that  there  exists  a  unique  solution  u  for  ( 1 )-  (2)  and  K  is  the  finite 
element  solution  with  linear  finite  element  spaces.  Then 


«  -  ||®iu,-p+t,i  +  Chm„.-r.2 


K 


2tn  i 


where  h  is  the  mesh  size  of  the  triangulations,  N  is  the  degree  of  polynomial,  K  is  the 
dimension  of  Gaussian  in  the  polynomial  chaos  expansion  and  C  is  a  constant  independent 
of  h,  N  and  K. 


In  many  cases,  the  wavenumber  is  treated  as  a  constant  and  a  random  noise  forcing  term  is 
added  in  Equation  (1)  instead  to  count  for  model  uncertainty.  We  have  studied  the  finite 
element  method  for  such  a  problem  and  obtained  rigorous  error  estimates(see  [4]  for  details). 


Conclusion  and  future  research 

We  posed  the  problem  of  optimal  control  and  design  optimization  of  acoustic  liner  under 
uncertainty  to  reduce  the  radiated  engine  noise  given  the  source.  Our  numerical  experiments 
show  that  the  choice  of  optimal  impedance  factor  results  in  reducing  the  far-held  noise  level 
significantly.  We  have  also  constructed  numerical  algorithms  to  solve  the  underlying  sto¬ 
chastic  partial  differential  equations  and  obtained  rigorous  error  estimates 


In  future  research,  we  plan  to  study  the  shape  design  of  the  fan  inlet  to  reduce  the  noise 
radiation  with  uncertain  wave  numbers.  We  also  plan  to  study  the  stochastic  problem  with 
the  wave  number  as  a  random  held,  not  just  a  random  variable. 
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